Dynamical Quark Effects on Light Quark Masses 
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We present results for light quark masses from a systematic lattice QCD study with two degen- 
erate flavors of dynamical quarks. Simulations are made with a renormalization-group improved 
gauge action and a mean-field improved clover quark action for sea quark masses correspond- 
ing to mps/mv ~ 0.8-0.6 and the lattice spacing a — 0.22-0.11 fm. In the continuum limit 



we find m*^*(2GeV) = 3.44 



+0.14 



MeV using the vr and p meson masses as physical input, and 



'(2GeV) 



MeV or 901^^ MeV with the K ot (f> meson mass as additional input. The 



quoted errors represent statistical and systematic combined, the latter including those from contin- 
uum and chiral extrapolations, and from renormalization factors. Compared to quenched results, 
two flavors of dynamical quarks reduce quark masses by about 25%. 



PACS numbers; 12.38.Gc, 12.15.Ff, 14.65. Bt 

Masses of light quarks belong to the most fundamental 
parameters of the Standard Model [1] , and yet their pre- 
cise values have been difficult to determine due to quark 
confinement. Lattice QCD provides a fundamental ap- 
proach to overcome this problem [2] since it enables first 
principle calculations of hadron masses as functions of 
quark masses. This approach has progressed consider- 
ably recently [3-5] through high statistics calculations 
that have allowed the continuum limit to be taken for 
quark masses, and the development of non-perturbative 
renormalization techniques for a reliable conversion of 
lattice bare quark masses to those in the continuum. 

These studies, however, have been carried out within 
the quenched approximation which ignores effects of sea 
quarks. A limitation, shown in Ref. [4], is that the 
strange quark mass cannot be consistently determined, 
with the values differing by 20% depending on the choice 
of meson mass taken for input. It has been suspected 
[2,3], furthermore, that dynamical sea quark effects siz- 
ably reduce the values of light quark masses. 

Clearly, systematic full QCD studies incorporating dy- 
namical sea quark effects are needed for progress in the 
determination of quark masses. A recent attempt has 
been reported in [6]. In this Letter we present results of 
our investigation in which the u and d quarks, assumed 
degenerate, are simulated dynamically while the s quark 
is treated in the quenched approximation [7] . 

Full QCD simulations are computationally much more 
demanding than quenched simulations. This problem can 
be significantly eased by the use of improved actions. 
Because of reduced cutoff errors, they should allow con- 
tinuum extrapolations from coarser lattices, and hence 
require smaller lattice sizes, and smaller computational 
costs, for simulations with the same physical lattice size. 



We employ improved actions both for the giuon part 
and the quark part. The gluon action consists of 1 x 1 and 
1x2 Wilson loops whose coefficients are determined by 
an approximate renormalization-group analysis [8]. For 
the quark part we choose the "clover" improvement of 
the Wilson action [9], adopting, for the clover coefficient, 
a mean- field value csw — P^"^^^- We substitute the one- 
loop result P = 1 - 0.140252 [8] for the plaquette P, 
which agrees within 8% with the values measured in our 
runs. The one-loop result of csw [10] is found to be close 
to our choice. 

The improved action described here was tested in our 
preparatory full QCD study [11]. We found that scal- 
ing violation in hadron masses is small with this action 
already at « 1 GeV, as compared to a^^>2 GeV 
needed for the standard plaquette and Wilson quark ac- 
tions. We therefore aim at a continuum extrapolation 
from simulations made at « 1-2 GeV. 

We make runs at three values of the coupling f3 = 
= 1.8,1.95,2.1 to cover this range, employing lat- 
tices of a similar physical spatial size La « 2.5 fm, as 
listed in Table I. For each (3, gauge configurations are 
generated by the hybrid Monte Carlo algorithm at four 
values of the sea quark hopping parameter Kg^a corre- 
sponding to the pseudoscalar (PS) to vector (V) meson 
mass ratio of Mps/My « 0.8, 0.75, 0.7 and 0.6. For 
each sea quark mass, we calculate hadron masses at five 
values of the valence quark hopping parameter Kyai cor- 
responding to Mps/Mv w 0.8, 0.75, 0.7, 0.6 and 0.5, 
taking unequal as well as equal quark mass cases. Masses 
are extracted from hadron propagators with the standard 
correlated fit- Errors are estimated with the jackknife 
procedure with a bin size of 50 trajectories, derived from 
an autocorrelation study. 
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For the Wilson-type quark action including the clover 
action, different definitions of quark masses lead to re- 
sults that differ at finite lattice spacing due to explicit 
breaking of chiral symmetry. Wc employ three defini- 
tions in the present work, checking consistency among 
them for reliability of results: (i) using the axial vector 
Ward identity, we define, 



lim 



2(p(i)p(o)) ■ 



(1) 



with P the pseudoscalar density and = A4-I-CAV4P 
the axial vector current improved to 0{a): (ii) an- 
other possibility, suggested by the vector Ward identity 
and naturally appearing in perturbative analyses, reads 
rn^^^a = — 1/kc)/2 where Kc represents the critical 
hopping parameter at which the PS meson mass van- 
ishes Mps{Kvai = Ksea = Kc) = 0; (hi) & third pos- 
sibility, suggested in [12] and denoted by rn^^^'^'^, re- 
places Kc in (ii) by the "partially quenched" critical value 
AjPQ where the PS meson mass vanishes as a function of 
the valence hopping parameter K^ai when Hg^a for sea 
quark is fixed to the physical point of u and d quark, 

MpsiKyal = ^; Ksea = Kud) = 0. 

We express the PS meson mass Mpg in terms of quark 
masses by a general quadratic ansatz of the form, 

MpgO? = bsTUseaa + byTTlyaia + Cs{mseaaf 

+ Cy{niyaiaf + Csvmseaamyaia- (2) 

Here ruyai = {mvai,! +m^ai,2)l'^ with myai,i{i = 1,2) 
the bare mass of the valence quark and antiquark of the 
PS meson, and msea the mass of sea quark. We mention 
that details vary depending on the definitions, e.g., terms 
depending only on nisea are absent for the case of AWI 
mass since with (1) Mps is zero for vanishing m^aU and 
a cross term oc myai,iO'mvai,20. is found to be necessary 
in the case of VWI mass. 

The vector meson mass My is written in a similar man- 
ner, adopting, however, PS meson masses as independent 
variables. We fit data with the formula 

Mya = A + BslJLseaa + B 



seo" 



+ ^vfJ'lal'^'^ + CsDMseaOMfaifl- (3) 



Here ^vai = (a*i + A*2)/2 represents the average of PS 
meson mass squared Hi = Mpg{Kvai,i, Kvai,u Ksea) made of 
a degenerate quark- antiquark pair of the valence hopping 

parameter Kyal,ii and j^sea ~ ^psi^seai Ksea] Ksea)' 

Fitting our hadron mass data with (2) and (3) we find 
reasonable results with x^/A'df in the range 0.6 2.3 (ex- 
cept for (2) for VWI quark mass at ^ = 1.8 for which 
X^/Nbf = 4.0). We then determine the bare lattice 
value of the average u and d quark mass rriudO. by fix- 
ing the valence and sea quark masses to be degenerate 
in (2) and (3), and requiring the experimental value for 



the ratio M^^/Mp = 0.1757. For the s quark mass TO^a 
we use either Mx/Mp = 0.6477 or M^/Mp = 1.3267 
while keeping the sea quark mass nisea at the value m„d 
determined above. The lattice scale is set using 
Mp = 0.7684 GeV as input. 

We convert bare quark masses calculated above to 
renormalized quark masses in the MS scheme at /i = 



1/a through 



vwi^ 



Zm{g )m^ 



vwifl and m^wi. 



{ZAig^)/Zp{g^))mf^^a. In these relations, the 0{a) 
improvement terms with the coefficients bm, bA and bp 
are also included. For renormalization factors and im- 
provement coefficients, including that for ca, one- loop 
perturbative values for massless quark [13] are used. 
For the coupling constant we adopt a mean-field im- 
proved value in the MS scheme appropriate for the RG- 
improvcd gluon action (see Table I for numerical values): 
gj^{l/a) = (3.648^^1x1 - 2.648H^ix2)/3/6 - 0.1006 -t- 
0.03149A^/ where measured values extrapolated to zero 
sea quark mass are substituted for the Wilson loops. 
The results for quark masses are run from n = 1/a to 
H = 2 GeV using the 3-loop beta function for Nf = 2 
[14]. Numerical values of quark masses at each /3 are 
given in Table II. 

In Fig. 1 we show our two-flavor full QCD results 



for 



(2 GeV) with filled symbols. The values for 



the three definitions, while differing sizably at finite lat- 
tice spacings [6] , tend to converge toward the continuum 
limit. A similar trend has been seen in the quenched data 
for the standard action [4] , reproduced in Fig. 1 with thin 
open symbols. 

For our choice of the improved action, the scaling vio- 
lation starts at 0{gj^{^)a) for the quark masses at the 
scale /U = 2 GeV. We tnerefore make a continuum extrap- 
olation linear in a. The results are given in Table II. The 
masses for the three definitions are consistent with each 
other at two-sigma level of statistics. Hence we carry 
out a combined linear fit, as shown in Fig. 1 by dashed 
lines, obtaining niud = 3.44(9) with x^/^of = 2.9 where 
the error is only statistical. The systematic error of the 
continuum extrapolation is estimated from the spread of 
values obtained by separate fits of data for the three def- 
initions. The fractional error thus calculated is given in 
Table III. 

This table lists our estimate for two more systematic 
errors that we need to incorporate. One is an uncer- 
tainty due to chiral extrapolations. We estimate this 
error from the change of the combined linear fit in the 
continuum limit when the quadratic term fj? in the vec- 
tor mass formula (3) is replaced by /i^/^ or cubic terms 
are included in the PS mass formula (2). Another is 
the error due to the use of one-loop perturbative values 
for the renormalization factors. As non-pcrturbative val- 
ues are not yet available, we estimate the effect of higher 
order contributions by recalculating masses while either 
shifting the matching scale from ^ = 1/a to fi = n/a 
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or using an alternative definition of coupling given by 
g^{l/a) = H^ix 1/3/6 + 0.2402 + 0.03149iV/ using only 
tne plaquette. 

Combining the statistical error and the systematic er- 
rors listed in Table III by quadrature to obtain the total 
error, we find for our final value, 

m^(2 GeV) = iAA+_^^\\ MeV (TV/ = 2). (4) 

Our full QCD result for the average u-d quark mass is 
considerably lower than our previous quenched result for 
the standard action given in [4] as 

m^(2GeV) = 4.57(18) MeV [Nf = 0, standard), (5) 

where the error is only statistical. In order to confirm 
that the decrease is a dynamical sea quark eff'ect, we 
carry out a quenched simulation using the same improved 
gluon and quark actions as for the full QCD runs. This 
simulation is made at 10 values of /3 chosen so that the 
string tension matches that of two-flavor full QCD for 
each simulated value of sea quark mass and for the chiral 
limit at /?= 1.95 and 2.1. 

Analyses leading from hadron masses to quark masses 
parallel those for full QCD. In particular we employ poly- 
nomial chiral expansions of the form (2)-(3), except that 
terms referring to sea quark masses arc dropped. As a 
cross-check we also make an analysis parallel to the one in 
[4], employing quenched chiral perturbation theory for- 
mulae, and obtain consistent results. 

We plot results of this quenched analysis with thick 
open symbols in Fig. 1. Good consistency is observed 
between the continuum values for the standard and im- 
proved actions. We also note that scaling violations are 
visibly reduced for the latter. Making a combined linear 
extrapolation we obtain in the continuum limit, 

m^(2 GeV) = 4.36tg;l^ MeV {Nf = 0, improved), (6) 

where the error is estimated in a similar way as for full 

QCD. From this analysis we conclude that the effect of 
two dynamical quarks is to decrease m^'^(2GeV) by 
about 25%. 

In Fig. 2 we show results for the strange quark mass 

m^^{2Ge\) determined from the K meson mass. A 
parallel figure obtained with the meson mass is given 
in Fig. 3. Using K* instead of cj) gives the same results 
within 1%. The strange quark is heavy enough so that 
the difference between Kg and k^'^ has only small effects 
on the VWI masses in full QCD. Employing combined 
linear continuum extrapolations in a (with x^/^DF = 1-3 
and 3.0, respectively), and estimating the error in the 
same way as for m„rf (see Table III for details), we obtain 

mf^(2 GeV) = 88tg MeV Mk input {Nf = 2) (7) 
= QOtfi MeV input {Nf = 2). (8) 



With (4) this gives mf^/m^f = 26(2) to be compared 
with 24.4(1.5) [15] computed from chiral perturbation 
theory to one loop. 

Similar analyses for quenched QCD lead to the results 

Mk input: 

mf^(2 GeV) = llOt^ MeV {Nf = 0, improved) (9) 

= 116(3) MeV {Nf = 0, standard) (10) 

input: 

mf^(2 GeV) = 132^^ MeV {Nf = 0, improved) (11) 

= 144(6) MeV {Nf = 0, standard), (12) 

where the values for the standard action [4] are also 
quoted for comparison. 

The quenched values for the standard and improved ac- 
tions are mutually consistent for each choice of the input. 
This confirms the existence of a systematic uncertainty 
of 20-30% in the value of in quenched QCD [4]. 

One of our important results for the strange quark 
mass is that this uncertainty disappears within an error 
of 10% by the inclusion of two flavors of sea quarks. The 
consistency reflects a closer agreement of the K — K* and 
K — (j) mass splittings with experiment in our two-flavor 
QCD results compared to the quenched case [7] . 

Another important result is that dynamical quark ef- 
fects reduce the value of m*^'^ significantly, from the 
range 110-140 MeV in quenched QCD to 90 MeV for 
two-flavor full QCD. It will be interesting to see whether 
the inclusion of dynamical effects of the strange quark it- 
self would decrease the value of mf'^ even further. This 
is an important issue to settle as our two-flavor results 
are already close to the lower bounds estimated from the 
positivity of spectral fimctions [16]. 

Clearly, establishing the values of light quark masses 
incorporating three flavors of dynamical quarks will be 
one of the main tasks of future lattice QCD calculations. 
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FIG. 2. Continuum extrapolation of the strange quark 
mass rus using Mk as input. Symbols have the same meaning 
as in Fig. 1. 
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FIG. 3. Continuum extrapolation of the strange quark 
mass TUs using as input. Symbols have the same meaning 
as in Fig. 1. 



a [GeV ] 



FIG. 1. Continuum extrapolation of the average up and 
down quark mass rriud for full QCD (filled symbols) and 
quenched QCD (thick open symbols) obtained with the im- 
proved action, and quenched results with the standard action 
(thin open symbols). Lines show combined fits linear in a. 
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TABLE I. Run parameters of two-flavor full QCD simulations. The scale a is set by Mp = 768.4 MeV. 
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X T 


Csw 


a(fm) 


La(fm) 




mps/mv for sea quarks : #traj. 






1.80 
1.95 
2.10 


3.168 
2.816 
2.565 


12^x24 
16^x32 
24^x48 


1.60 
1.53 
1.47 


0.215(2) 
0.155(2) 
0.108(1) 


2.58(3) 
2.48(3) 
2.58(3) 


0.807(1): 6250 
0.804(1): 7000 
0.806(1): 4000 


0.753(1): 5000 0.694(2) : 7000 
0.752(1): 7000 0.690(1) : 7000 
0.755(2) :4000 0.691(3) :4000 


0.547(4) 
0.582(3) 
0.576(3) 


5250 
5000 
4000 



TABLE IL Renormalized quark masses (in MeV) in the MS-scheme at jj, = 2 GeV for each (3, and in the continuum obtained 
by linear fits in a (for these fits x^/^df is also given). All errors are statistical. 



/3 


'"-ud 


VWI.PQ 






^ywi,PQ 


nis (a) 


my^^ (0) 


my^^'P-^ (0) 




1.8 


2.277(27) 


4.183(42) 


3.322(37) 


102.92(92) 


104.54(93) 


88.0(1.0) 


129.1(2.2) 


130.7(2.2) 


113.9(2.4) 


1.95 


2.489(38) 


4.064(43) 


3.321(38) 


100.65(98) 


102.08(99) 


87.2(1.0) 


123.1(1.7) 


124.5(1.7) 


109.8(1.7) 


2.1 


2.966(55) 


3.816(47) 


3.344(46) 


95.6(1.1) 


96.4(1.1) 


87.0(1.2) 


108.0(2.2) 


108.8(2.2) 


100.0(2.2) 


a 


3.47(10) 


3.50(10) 


3.36(9) 


89.4(2.3) 


89.5(2.3) 


85.8(2.4) 


90.1(4.9) 


90.3(4.9) 


88.1(4.9) 


XV^DF 


10.8 


2.4 


0.07 


2.1 


2.7 


0.03 


6.0 


6.5 


2.4 



TABLE in. Contributions to total error in continuum limit. 



statistical chiral Z-factor cont. cxtrap. 



mud 


+2.6% 


+ 1.2%o 


+2.3%o 


+1.7% 




-2.6% 


-2.3% 


-5.0% 


-2.3% 


rus (/f -input) 


-h2.4% 


-1-1.6% 


-h2.2% 


-hl.4% 




-2.4% 


-2.2% 


-5.6% 


-2.8% 


ms (0-input) 


+4.8% 


+1.5% 


+1.7%o 


+0.9% 




-4.8% 


-7.6% 


-6.9% 


-1.6% 
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